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It is shown that coherent scatterings by an ordered vortex lattice are critically enhanced for 
quasi particles moving in cyclotron orbits on the Fermi surface through vortex core regions, thus 
generating significant quasi-periodic oscillating contributions to the SC free energy as a function of 
the inverse magnetic field. The mean frequency of the oscillation provides a fingerprint of the vortex 
lattice geometry. Vortex-lattice disorder, tends to suppress this oscillatory component. 
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Many reports on observation of de Haas van Alphcn 
(dHvA) or Shubnikov-de Haas (SdH) oscillations in the 
mixed (vortex) states of strongly type-II superconduc- 
tors, have shown additional damping of the oscillations' 
amplitude in the vortex state with respect to the normal 
state signal (see, e.g*^). In all these experiments the 
appearance of the extra damping has been found struc- 
tureless, reflecting only some average effect of the vor- 
tex matter on the fermionic quasi particles. Extracting 
valuable information about the vortex state from these 
data requires development of a comprehensive quanti- 
tative theory of pure strongly type-II superconductors 
at high magnetic fields and low temperatures, a the- 
ory which docs not exist today even within the conven- 
tional BCS framework^. Recently, using a theoretical ap- 
proach based on an exact perturbative Gorkov-Ginzburg- 
Landau theory^, it was found, however, that on certain 
cyclotron orbits fermionic quasi particles (QPs) are sin- 
gularly coupled to, and coherently scattered by the vortex 
lattice, resulting in a new type of magnctoquantum os- 
cillations superimposed on the usual dHvA oscillations. 
In the present paper we analyze characteristic features of 
the discovered oscillatory effect and discuss their experi- 
mental feasibility. The kinematical condition controlling 
the predicted effect is shown, by means of a simple model, 
to be operative beyond the limitations of the perturba- 
tion theory. 

Following the theory presented in Ref£ we consider a 
2D strongly type-II (neglecting the effect of SC screen- 
ing currents) superconductor in a perpendicular uniform 
magnetic field H =Hz. Generalization to isotropic 3D 
systems is rather straightforward. It is assumed that 
the superconductor can be described by means of BCS- 



Hamiltonian for the usual singlet s-wave electron pair- 
ing (electron spin is neglected for the sake of simplicity) . 
Within mean-field approximation (but not in a fully self- 
consistent manner) the order parameter is described by 

a general vortex lattice state, A(r) = I 2? J A <po( r )) 
written in terms of a discrete set of ground-state Landau 
orbitals (in symmetric gauge): 

<po(x,y) = e ix y J2 n e- %9n2+l ' 1 " x -(y +q "/ 2 '> 2 , where q n = 
—n = qc,n , n = 0,±1,±2, ... with the lattice spac- 

ing a x along the x - axis and a 2 , = tt/ \J 1 — (6/ir) 2 . 
For the Abrikosov triangular and square lattices: 9 = 
7r/2,and 9 = respectively. We use dimensionless space 
coordinates measured in units of the magnetic length, 
an = v ch/eH. The amplitude of the order parameter, 
A 2 , = S- 1 J d 2 Y t |A(r,)| 2 , where S = ttN and N is the 
number of vortices, is treated as a variational parameter 
with respect to the SC thermodynamic potential (TP) 
Q sc (Ao), which can be written as a Taylor expansion in 
A : 

tt sc (A ) = sf^) + vH2!f> 2n (A ) (1) 



where n 2n (A ) = Nk B T I 2n , g BCS is the BCS 

coupling constant and uj c = eH/mc. The quartic term, 
I4 , is the leading contribution to the SC free energy 
influenced by the vortex distribution^. An exact ex- 
pression derived ir£ for a given Matzubara frequency 
uj v = irksT (2v + 1) /h, reads: 



dT 1 dT 2 dT 3 dT 4 e"^" (Tl+T2+T3+r - l) - m ^ Tl - T2+T:i - Ti) ^ (2) 

a\ + 012 + Q'3 + «4 



where vj v = lu,j/uj c and rip = Ep/huj c — 1/2. Here j3 (7) = ^2 g Pg (7) is a function of the 4-electron variable 
7 = ot_i = 1 - e*^, e_i = (-1V+ 1 , and: 
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FIG. 1: (color online): Left panel: The "erratic" function w (uf) (solid blue line) plotted together with the characteristic 
smooth function S (np) (dashed red line) in the interval 20 < uf < 100. Right panel: Number of maxima, iV max {y/npj, of the 
characteristic function H (jif) for a square vortex lattice (red line) and for a hexagonal vortex lattice (blue line). 
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where G is a reciprocal vortex lattice vector. The structure function /3 (7) controls the coupling between the 
two pairs of electrons involved and the vortex lattice, where 7 is a periodic function of the individual electronic 
"time" variables Tj, which reflects the underlying electron cyclotron motions. Note that /3 (7 = 0) is the well known 
Abrikosov parameter /3a for an arbitrary vortex lattice geometry^. The most remarkable feature of j3 (7) seen in 
Eq.Q is associated with the dual singular points at 7 —> ±1, where Tj, approach, respectively: T1 = t 3 — > 0,t 2 — s- 
nn — t,T4 — > nir + r , or: n — > mr — r, t 3 — > nir + r, r 2 = t± — > , with r being an arbitrary real number in 
the interval: — it < t < it, and n = 1,2,.... The electrons at such highly correlated pairs of cyclotron orbits are 
singularly scattered by the vortex lattice, yielding only purely harmonic contributions to the SC free energy in the 



dHvA frequency F = ufH since under these conditions: 



if(ti-t 2 +t 3 -t 4 ) 
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The pair of electrons 



whose positions on the singular cyclotron orbit coincide (i.e. (1,3) at 7 — > 1 and (2,4) at 7 — » —1 ) undergo local 
mutual scattering and so exchange many G- vectors through the vortex lattice during the scattering process, while 
those electrons moving coherently on a large cyclotron orbit in opposite directions ( i.e. (1,3) at 7 — > —1 and (2,4) 
at 7 — ¥ 1) are mutually scattered over a long range of the vortex lattice, and so do not exchange momentum. Thus, 

the leading purely harmonic contributions to j3 (7) are: /3g=o (7) = § (177 + 1^7) ' anc ^ / ^ G W d 2 G, which yields: 

h rr-^ — for 7 — > 1 , and 7^— for 7 — > —1. 

2 l— 7 ' ' 1+7 ' 

Significant deviations from the harmonic part are developed slightly away from the singular points. They can be 
estimated by expanding (1 — 7) / (1 + 7) about the singular points, e.g. for 7 = 1: (1 — 7) / (1 + 7) ~ — + 

+ £3) A 6 > with £l = ( T i +t 2 +t 3 + r 4 ) /4 - it/2 , I2 = (ti - r 2 + r 3 - r 4 ) /2 + tt , £ 3 = n - t 3 , and carrying 
out the Tj-integrations. Focusing, for simplicity, on the deviations from the first harmonic it is found that: I4 — > 
Ree 2 OT in 1 P e -2^„ j°° d^e-^-^h (it) , where: 
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It is clear that the dominant contributions to the in- 



tegral in Eq.@ originate from reciprocal lattice vectors 
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FIG. 2: (color online): A square lattice of vortex cores (small red circles), shown sweeping across a pair of classical orbits 
occupied by a quasi particle-quasi hole pair (upper and lower solid curves respectively) in a magnetic field whose magnitude 
increases successively from the upper left panel through the lower right panel. Note that in the ground state the occupied 
(low-energy branch) cyclotron orbit is locally distorted paramagnetically in passing through vortex core regions, wherease both 
orbital branches of a quasi particle-hole excitation are distorted diamagnetically. 



satisfying: |G| ps 2\J2np, namely having length close 
to the Fermi surface diameter. Furthermore, due to the 
large values of np and the discrete nature of G (which 
are measured in units of a^ 1 ), the integration over £2 
yields "erratically" oscillating, quasi-periodic function of 
y/nF, superimposed on the usual dHvA oscillation, ex- 
pressed by the factor e 2mnp of I4 in our simple spheri- 
cal Fermi surface model. As will be illustrated below, 
the corresponding oscillating envelop function contains 
vortex-lattice structural information which can be ex- 
tracted from experiment by an appropriate data analysis 
(see Fig.l). 

The final result for the first dHvA harmonic of the SC 
TP, up to fourth order, can be written in the form: 
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where fin' 1 " 1 is the corresponding normal state free energy 
extrapolated to the vortex state, and w(jif), shown in 
Fig.l, represents the "erratically" oscillating component. 

The influence of vortex-lattice disorder on the SC 
free energy in the white noise limit can be determined 
by invoking the general expansion of the state function 
tpo{%> y) in terms of Landau orbitals wave functions with 
randomly selected coefficients. Averaging (3 (7) over real- 
izations of these coefficients it can be easily shown that: 
(j3 (7)) — > + yA^j which is just the purely harmonic 
component of f3 (7). In this limiting case, only incoher- 
ent scattering processes by the vortex matter contribute 
to the SC TP, and the final result, up to fourth order, 
is obtained from Eq.([5]) by taking w (np) — > 0, i.e. very 



close to the well known Maki expression^, as expanded 
to the same power in Aq. 

It is interesting to note that the kinematical condition, 
|G| = lyjlrip , is equivalent to the real space condition, 
|R| = 2y/2np (where R is a vortex-lattice vector mea- 
sured in units a# ), which is just the condition for the cy- 
clotron orbit at the Fermi energy to pass through a vortex 
core. The "erratic" oscillations are closely related to such 
cyclotron orbits since the latter are strongly distorted 
during their passage near a vortex corei The influence 
of a "sweeping" vortex lattice (generated by a sweeping 
magnetic field) crossing a quasi-particle quasi-hole pair 
of cyclotron orbits may be studied qualitatively by con- 
sidering a highly simplified model of a classical charged 
particle moving in two dimensions under a perpendicular 
magnetic field H = H'z and in the presence of a "pair- 
■ potential", |A(r)| , with the Bogoliubov dispersion rela- 

h 2 k 2 



tion e (k) = ±y/t* + |A (r)| , where & = ^ - cf*. 
An important advantage of this model, despite its obvi- 
ous shortcomings, is that it can be implemented without 
invoking perturbation theory with respect to A (r) . The 
corresponding coupled Lorentz-Bloch velocity equations: 
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can be easily integrated to yield the trajectory equation: 



(kp/hY 



where h = H / ' H C 2 and ro is an integral of motion de- 
termined by initial conditions. Obviously, there are two 
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types of solutions for energies above and below the Fermi 
energy, corresponding to the Bogoliubove particle-hole 
pair. Examples of classically allowed, distorted cyclotron 
trajectories, which correspond to real solutions r ((/>) , 
are depicted in Fig. 2. They are restricted to QP ener- 
gies above the threshold Tq = 2 max | A (r, (f>)\ /ftw^ti 1 . 
Below this threshold there are complex solutions corre- 
sponding to tunneling of QP orbits between vortex cores 
(not shown). For increasing magnetic field strength h at 
a fixed QP energy, the lattice of vortex cores, shown in 
Fig. 2 successively from the upper-left panel through the 
lower-right panel, sweeps across the QP cyclotron orbits, 
which undergo a strong local distortion upon crossing a 
vortex core. The appearances of such local distortions are 
closely related to the "erratic" parmagnetic-diamagnetic 
oscillations of the function w (jif) shown in Fig.l. A 
characteristic feature of w (njr) to be exploited for distin- 
guishing between different vortex-lattice structures can 
be conveniently defined in terms of the smooth function 



S(tif) = aV^exp 



-i|G| 2 (l-8n F /|G| 



— b (see 



Fig.l, left panel), derived with the help of Eq.@ by per- 
forming the integral J °° d^ie^ 4 ^"^ 1 13 (^ij in the station- 
ary phase approximation. Here a and b are adjustable 
parameters. The number of maxima, iV max (y/nF~) of 
S(n^) appearing bellow a given value of np offers the 
desired characteristic feature (sec Fig.l, right panel). Its 
slop is a measure of the rate at which vortex cores in 
Fig. 2 cross the cyclotron orbit at the Fermi energy, and 
so uniquely characterizes the point symmetry of the un- 
derlying vortex lattice. 
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